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ABSTRACT 


This final report describes the aims of this research and the experimental 
equipment set up at the Research Annex of Texas A&M University. A review 
of theoretical considerations is presented together with a summary of the attempts 
made to obtain a numerical solution and the difficulties encountered. 

The final section of this report describes new work concerned with the 
experimental determination of the chordwise pressure distribution during 
operation at a tip speed of 500 ft /sec. Chordwise integration of the pressures 
leads to the spanwise load distribution and further integration would give the 


axial thrust. 
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INTRODUCTION 

i 

The work reported here has been carried out under NASA Grant NSG 669 
initiated under the direction of Messrs. A. E. Cronk and J. C. Brusse. In 
submitting the renewal proposal in February 1967 , the investigators were 
Messrs. A. E« Cronk, C. F. Kettleborough and J. C. Brusse. At the end of 
the year Prof. A. E. Cronk was no longer a member of the team. 

The problem to be studied was that of improving the performance estimation 
of propellers operating in the heavily loaded static thrust condition. The 
design of aircraft propellers by analytical methods, using the classical 
propeller theory of Goldstein has not been completely successful, particularly 
when performance optimization for the static thrust condition is attempted. 

The origin of the deficiencies in the analytical method is not apparent and 
has caused considerable concern in the development of propellor driven V/STOL 
aircraft. 


EXPERIMENTAL CONSIDERATIONS 

In order to assess the Goldstein Theory as it applies to propellers oper- 
ating in the static thrust condition, it was considered necessary to in- 
augurate an experimental investigation to evaluate the performance of propellers 
alone and then to attempt correlation of analytically developed performance 
data with the experimentally generated data. It was also considered necessary 
to acquire a more complete understanding of the actual airflow phenomena in the 
vicinity of the propeller in order that the basic assumptions involved in the 
existing analytical methods can be assessed. Critical evaluation of the several 
propeller performance analytical methods Indicated that a better knowledge of the 
radial loading distribution as well as the inflow of air into the propeller was 
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required. Consequently it appeared necessary that the experimental apparatus 
permit measurement of the air-flow at the propeller plane of rotation and 
perhaps the actual measurement of the radial load distribution on the blade 
itself. It was felt that if the flow patterns in the real fluid are understood 
more completely, the development of an analytical method would be much more 
easily attained. The ultimate aim of the research was the development of an 
analytical method which will yield performance estimates that more nearly 
compare to the actual performance observed. 

The program called for the construction of a propeller static thrust 
dynamometer which would be adequate for conducting propeller performance ex- 
periments as well as the airflow studies. In addition to the static thrust 
dynamometer, a program of study and development of analytical methods useful in 
the prediction of propeller performance was to be carried out. 

A propeller static thrust: dynamometer designed to minimize wake and in- 
flow interference has been constructed in a large airplane hangar at the Texas 
A&M Research Annex. The size and configuration of the dynamometer was deter- 
mined after careful consideration of the following desirable features: 

1. Test propellers must be isolated to the greatest degree possible. 

2. The dynamometer must be capable of operating model propellers large 
enough to permit evaluation of Reynolds' Number effect. 

3. The facility should be installed inside a building to avoid the 
vagaries of the weather. 

With these factors in mind a configuration was conceived wherein the propeller 
is mounted so that it rotates in a horizontal plane with the line of thrust 
acting vertically upward. The propeller is rotated by a vertical shaft which 
turns inside a 7.5 inch diameter tubular steel housing. The housing has four 
slanted legs welded to its lower end which are bolted to specially constructed 
reinforced concrete piers eleven feet deep in the earth. The propeller plane 

ORIGINAL PAGB m 
OF' POOR QUALITY 



3 . 


of rotation is 20 feet above the floor of the hangar and 10 feet below the roof 
joists. The curved hangar roof provides about 15 feet of additional over- 
head space that contains the roof supporting trusses. The nearest walls are 
in excess of 40 feet from the propeller. While the desirable condition of 
testing in an infinitely large air space has not been totally realized, it is 
believed the installation approximates this condition within practical limits. 

Preliminary calculations indicated that approximately 400 horse-power would 
be required to operate 5-foot diameter, 4-way, high activity factor propellers 
at tip speeds of 900 feet per second. In addition to providing these values 
of power and speed, the power source must be easily controllable through a wide 
speed range. In the interest of economy, internal combustion engines were 
chosen for the power source instead of an electric motor drive. Two 275-HP 
Chris Craft V-8 Marine engines are installed to drive the vertical shaft 
through a differential gear box. Couplings are installed at each engine out- 
put shaft to allow for installation misalignment and vibration isolation. An 
overrunning clutch is mounted between the differential gear box output shaft 
and the main drive shaft to provide protection for the propeller balance in 
the event of mechanical failure in either of the power trains or the gear box. 

A splined axial-slip coupling is provided between the upper end of the main 
drive shaft and NASA provided propeller shaft to allow for thermal expansion 
and contraction of the shafts. A NASA strain— gaged propeller balance is mounted 
upstream of the thrust bearing and a propeller hob adapter shaft and propeller 
are then mounted directly to the balance. The electric signals from the bal- 
ance are routed through a slip ring assembly and thence through conduit to 
amplifiers and the read-out gages. A 60- tooth steel gear is mounted, on the 
propeller shaft between the balance and the thrust bearing to provide im- 
pulses through a magnetic pickup to an electronic tachometer. A small control 
house was built to house the power controls and the data acquisition equipment. 



The slipring assembly originally used for transmission of the propeller 
balance strain gage signals caused considerable trouble in the initial stages 
of operation. Brush wear seemed excessive and because of inadequate air 
blast scavenging, the displaced brush particles were deposited around the 
inner walls of the brush holding assembly. These deposites provided variable 
resistance circuits between adjacent brushes and resulted in large fluctua- 
tions and zero shifts in the strain gage bridge output circuits. Cleaning 
of the slipring assembly was required after about 5 hours of operation and 
because of the necessity to recalibrate the data acquisition system after 
each disassembly the down-time was at least three days. 

A slipring assembly having adequate cooling and scavenging provisions 
was installed during the month of July 1965 and has operated satisfactorily 
and without attention since that time. Some zero-shift between the beginning 
and end of each test run continues to be manifested, especially during the 
first 10-15 minutes of operation while the slipring assembly is cold. It has 
been observed that if the dynamometer is run for 20 or 30 minutes prior to 
taking any data, the zero-shift is reduced to insignificant values. With 
careful calibration and adjustment of the bridge power supplies, It appears 
possible to generate data which is repeatable within 2 % of the measured per- 
formance of the propeller for C greater than .025. 

P 

A typical performance test for 5 foot diameter propellers requires the 
determination of thrust and torque for a range of rotational speeds from an 
arbitrary low value of 1500 RPM to the maximum allowable for the propeller be- 
ing tested in approximately 200 RPM increments. The blade is set at a pre- 
determined blade angle for each test run, the usual practice is to vary the 
blade angle at the reference station from zero to 25 degrees in 2.5 degree 
increments. Various performance parameters are calculated from the basic 
values of thrust, torque and rotational speed for each condition and blade 


angle tested. 
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In order to provide visual evaluation of the airflow surrounding the 
propeller, a smoke generator with appropriate nozzles has been developed. In 
addition to the smoke generator, a slavestrobe unit was purchased which permits 
a gas discharge flashtube to be synchronized at any desired propeller 
speed. The usual flow visualization technique required the light to flash 
once per revolution of the propeller. With the smoke generator discharging 
a filament of smoke at a prescribed point while being illuminated by the 
flashes of the slavestrobe, an observer is able to get an idea of the relative 
airflow velocity at various points as well as the direction of airflow through 
the propeller. In addition to visual observation of the airflow phenomenon, 
photographs with 3,000 ASA Polaroid film have been made. These photographs 
show the general inflow pattern of air into the propeller as well as the vortex 
system that is formed in the wake of the propeller. The helical path of 
the airflow in the propeller slip stream also is observable. Thus the intro- 
duction of smoke at various points around the rotating propeller has provided 
a three-dimensional evaluation of the flow phenomenon associated with the 
propeller operating in the highly loaded static thrust condition. 

THEORETICAL CONSIDERATIONS 

For many years the theory developed by Goldstein^ has provided a 
reliable means of calculating propeller performance. Modified by experience 
factors and empirical corrections, this theory has been used intensively in 
the design of propellers, largely for climb and cruise conditions of operations. 
Attempts to optimize propeller designs for static thrust operation have in- 
dicated discrepancies between theoretical performance and actual performance 
achieved when the Goldstein theory is applied. This research is directed to- 
ward investigating these discrepancies and attempting to develop improved 
methods of calculating propeller performance for static thrust condition. 
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Possible discrepancies between the mathematical model postulated in the 
Goldstein theory and the physical flow situation are apparent. 

1. Some of the power put into the propeller is lost in generating the 

trailing vortex field of the propeller wake. The theory assumes, 

( 2 ) 

following Betz , that the load distribution on the propeller blade 
is such as to minimize this power loss for a given thrust. In the 
heavily loaded static thrust condition there is no assurance that 
such an optimum load distribution is achieved. 

2. As stated by Goldstein, the theory is valid only for small values of 
the ratio of the displacement velocity of the assumed helical vortex 
sheet to the velocity of advance of the propeller. Since this ratio 
becomes infinite in the static thrust condition it appears to be 
overly optimistic to expect the theory to hold with good accuracy in 
this situation. 

3. The assumption of the propeller wake being in the form of helical 
vortex sheets of uniform pitch will be violated in the region immediately 
downstream of the propeller. In this region, herein called the near 
wake, the flow is accelerating in the axial direction and hence, the 
pitch is Increasing. The near wake will be dominant in determining 

the inflow velocity field associated with the trailing vorticity. In 
addition, the instability of the vortex sheet and the tendency to 
roll up into discrete vortices will affect the inflow. 

The interdependence of the blade load distribution, the trailing vortex 
field and the inflow velocity field poses a difficult problem that may not be 
amenable to a direct analytic solution. 

Two methods of calculating the inflow velocity fields were selected for 
use in this investigation. They are the method of Iwasaki (3) which is an ex- 
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tension of methods presented by Kucheman^ and the method of Both 

methods involve an initial knowledge of the blade loading and so lend them- 
selves to iterative procedures. In the experimental part of this program, the 
velocity fields were measured. 

In the approach used by Wu the propeller is regarded as an actuator disk, 
normal to a uniform steady flow of an incompressible, inviscid fluid, Re- 
ferred to a cylindrical coordinate system (r,9,Z) with the Z axis coincident 
with the axis of rotation, positive direction downstream, the fluid velocity 
components are (u,v,w) respectively in the radial, tangential and axial direc- 
tions. Since the flow is steady and axisymmetric a stream function can be 
expressed as a function of r and z from which two of the fluid velocity com- 
ponents u and w can be determined as 


u= 


_L 3^ 

r 1 


v\/ - - 4 -~ ^ 

r 2>r 


a) 


The tangential velocity component v, must be specified separately. Con- 
siderations of conservation of angular momentum of the fluid show that the 
product vr must be a function of Y in the slipstream and zero elsewhere. 
Writing the equations of motion of the fluid in terms of the stream 
function results in 

• in the slipstream 


vr 


0 elsewhere 


(2) 


^ _ i&P . (Jlr v -vr)ivr, (i?r-vr) r ^ r | vr (3) + l ^< ‘ vr 

^ r*r ‘ % V S . M 


where Jc is the constant angular velocity of the propeller. 
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This constitutes a set of non-linear, partial differential equations for 
q, and v. Two difficulties arise in obtaining solutions to these equations; 
namely, the nonlinearity apparent in the terms of the right-hand side of (3) 
and the fact that the slipstream region is not known. 

These equations can be written in dimensionless form normalizing on the 
radius R as the unit of length and the remote uniform velocity, W. Utilizing 
differential operators that apply in a linear fashion and a perturbation stream 
function, ^ , it is possible to express the stream function Ip in a non-linear 
integral equation with the domain of integration not defined. 

Since the domain of integration is determinable from the stream function 
a method of successive approximations is suggested. The problem of existence 
and uniqueness of solutions so obtained is mentioned by Wu but not treated 
by him. In addition, in the numerical calculation procedures some problems of 
numerical instabilities and lack of convergence in the interative procedures 
have appeared. 

The approach used in this study is to retain the equations developed in 
dimensionless form, normalizing on W, the remote velocity. Attempts were 
made to extrapolate to the static thrust condition of W=0, in order to facilitate 
the numerical calculation procedure which involves an initial assumption of a 
stream function corresponding to the free stream flow undisturbed by the 
actuator disk. For the static thrust condition this would reduce to the 
trivial case of W=0 at all points or tp = constant, which will not provide 
the required starting point for the calculation. 

A program for calculation of the velocity field was developed following 
exactly the calculations described by Wu. This is an iterative procedure. 

Wu states that the second iteration gives results that are in general suffi- 
ciently accurate. However, it was found that an instability developed in the 
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third iteration that casts doubts on the accuracy of the second iteration. 

It is therefore obvious that the second iteration does not produce a stable 
answer. To this date no numerical results have been published using Wu's 
method. Because of the lack of convegence of Wu's method of calculation an 
alternative approach was considered in which equation (3) was cast in the form. 


C - 



) 

op 


> ( - j 

i 


^ER^ciENTS: 

< 




^ J 1 

1 
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To specify the problem it is necessary to define the manner in which v 
varies in the slipstream. A simple form was originally defined as v=kr 
(where k is an experimental constant). The coordinates were made non-linear so 
that a large number of grid points were situated along the axis and near the 
propeller plane. 

Other forms of the induced tangential velocity v were considered; these 
elude: 

i n 

v » kr 



v » Ar a (1-r) 6 
v * A(r-r c ) n (l-r) m 


However, the shape of the slipstream downstream of the propeller is not 
known and no method was found that allowed the shape to be determined. It is 
known that the angular momentum in the slipstream must remain constant. 

Several models based on this principle of conservation of angular momentum 
were tested as part of the overall numerical solution. In all cases the shape 
of the slipstream varied each iteration and could not be stabilized. Another 
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method was to initially guess the shape of the slipstream and then check for 
consistency of angular momentum at each section. A correction to the 
slipstream diameter could then be estimated so that angular momentum would be 
conserved. However, this approach did not converge. Inspite of all efforts 
no method was found which satisfactorily calculated the slipstream shape down- 
stream and upstream of the propellor. 
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DETERMINATION OF BLADELOAD DISTRIBUTION 

Propeller Blades 

Two propeller test blades were made in accordance with the geometry 
described In figure 1. The unusual blade thickness was dictated partly 
by the requirement for pneumatic tubing to be inlaid into the blade surfaces 
and partly because bending and torsional stiffness was considered important. 

Since it was planned to conduct the tests at relatively low tip speed to stay 
away from compressibility effects, the 10% blade tip thickness was not considered 
detrimental. 

The blades were milled from solid 7075-T651 aluminum bar stock, both 
blades being copied from a master blade by use of a milling machine hydraulic tracing 
attachment. Construction of the master blade involved the development of a set 
of airfoil templates for each spanwise orifice station. Each template was 
then copied in the master blade stock by use of the milling machine tracer 
attachment. The designed twist distribution was machined into the master blade 
by rotating it about its spanwise axis (.50 chord) to the required angle for 
each station template. The spaces between the spanwise stations were rough 
milled for fast material removal and then hand filed to a smooth transition 
surface. The final master blade shape was achieved by alternatively filing 
and checking both chordwise and spanwise profiles with a recording profilometer 
until the finished blade was within .005 inches of the designed contour. 

After tracer milling two test blades from the master blade, spanwise 
grooves one-sixteenth wide were milled into both the upper and lower surfaces to 
provide space for inlaying the pneumatic tubing. Figure 2 presents the locations 
and the tube numbering for each blade. The grooving operation required the 
development and manufacture of a spacial cutter depth-limiting device and 
a holding fixture that provided blade rotation about its spanwise axis while 
the milling operation was in progress. 
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A paste type epoxy adhesive was used to bond one— sixteenth inch diameter 
hard brass tubing into the machined grooves. The paste was also bridged over 
the exposed tubing and provided a filler in the groove for subsequent filing 
to the blade contour. The shank end of each brass tube was fitted with a length 
of flexible plastic tubing which was then routed through appropriate holes in both 
upper and lower surfaces and out the end of the blade shank through a hole 
drilled along the blade spanwise centerline. After careful leak testing, all 
holes and pockets were expoxied over with more of the paste epoxy which, upon 
achieving complete curing, was also filed to the finished blade contour. 

All of theblade orifices were drilled into the brass tubing with a 0.015 
inch drill. Brass chips from the drilling operation were blown out the tip 
end of each tube with compressed air through the flexible tubing end. A piece 
of 1/32 inch diameter piano wire was used to ream each tube to remove any burrs 
at the junction of the drilled orifices and the tube bores. Great care 
was taken with the orifices to assure clean, round holes with as near identical 
lips as possible. 

After the blades were completed and polished, each one was weighed 
to determine its; weight and center of gravity. Only minor weight adjustments 
were required to bring both blades to near identical weight and c.g. location. 

Figure 3 presents a sketch of a test blade showing schematically the tubing 
layout and orifices. 

Propeller Hub 

The blade shanks were mounted in a ground adjustable, 4-way hub machined 
from AISI 4340 steel alloy, heat treated to 125,000 psi tensile strength. 

It was designed for flange mounting to the short shaft of the TAMU propeller 
static thrust dynamometer. A three- inch bore was provided along the axis of the 
hub to house a scanivalve and the flexible tubing from the blade orifices. A fiber- 
glass spinner was constructed to provide additional protected space for the hub- 
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mounted data acquisition electronic equipment. 

Scanivalve and Associated Equipment 

A model 48D3 Scanivalve with a specially modified stepping motor was used 
for scanning the blade orifices during the test runs. A mounting fixture was 
devised that positioned the Scanivalve axial centerline on the dynamometer 
shaft axis of rotation, resulting in the Scanivalve pressure transducer rotating 
about an axis normal to its diaphragm plane. Because of the symmetry of the 
transducer and the relatively slow rotational speed attained, the centrifugal 
loading in the transducer diaphragm plane was quite small. In use, there was 
no measurable error resulting from the rotating transducer. 

Special electronic devices were added to the rotating data acquisition 
equipment in the hub to minimized the effects of unsteady DC voltage transmission 
across the slipring assembly. Power for the Scanivalve was routed from a 
DC power supply in the control house to the stator side of the slipring 
assembly mounted on the dynamometer short shaft. Figure 4 presents a cut-away 
sketch of the propellor hub assembly. A miniaturized transducer excitation sub- 
regulator was mounted in the hub spinner along with the transducer balance and 
calibration device. The output from the transducer was changed to a FM signal 
in a miniaturized sub-carrier oscillator also mounted in the spinner. The 
FM data signal was transmitted through the sliprings and back Into the control 
house to a discriminator which changed it back into the original analog signal. 

A signal conditioner and amplifier were provided for the position encoding 
device associated with the Scanivalve. Figure 5 Is a block diagram of the data 
acquisition equipment. 

Data Acquisition and Redording 

The flexible plastic tubing extending into the hub bore from each propeller 
blade shank was connected to the Scanivalve tubulations in accordance with the 
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tubing numbering scheme shown in Figure 2. Scanivalve tubulations 41 through 
44 were used in the determination of the transducer reference pressure. 

In operation the Scanivalve was always started from the "home" position 
and stepped through each tubulation position until a complete scan of the pressure 
tubes was finished. For each Scanivalve position three separate pieces of data 
were recorded: propeller speed in rpm, Scanivalve position digits and the transducer 

output as a signed voltage. 

Propeller rpm was sensed by a Hewlet-Packard 55RA digital counter set to count 
for one second the pulses generated by a magnetic pickup operating near a 60-tooth 
gear that turned with the dynamometer shaft. The Scanivalve position was determined 
by a special counter/controller that counted the pulses furnished by the Scanivalve 
position encoder during each step. The pressure transducer discriminator output 
was measured by a HP2104C integrating digital voltmeter. All of these counting 
and measuring devices provided visual digital display and an encoded BCD 
output. 

The HP2515A high speed digital output scanner is a programmable switching 
device providing options for controlling up to nine data sources and a data 
recording device. As used for this test, the three data sources sensing rpm, 
position, and pressure were connected to the scanner with appropriate cabling 
while a HP 562A paper tape digital printer was used as the recording device. 

Testing 

A typical test began with checking or resetting the propeller blade 
angles. Each blade is set with an inclinometer at the same azimuth position 
to assure that each blade has the same setting. All pressure orifices are 
taped over with one mil mylar tape except the spanwise station on each blade to 
be measured. The transducer was electrically zeroed and the calibration resistance 



Inserted to check the bridge circuitry. A pressure scan was performed with 
the propeller stationary — all transducer values for all the orifices should 
be identical and zero. Any reading shown to be different usually indicated 
a scanivalve problem or an obstructed orifice or tube. 

After the "zero check" the dynamometer was started and brought to a speed 
between 1500 and 2000 rpm and allowed to run at that speed for 15 to 20 minutes 
to bring all mechanical and electronic equipment to normal operation 
temperature. Immediately after stopping the dynamometer, another "zero check" 
scan was performed to determine if any partial obstruction existed in any 
of the pneumatic equipment. If there was no indication of pneumatic obstruction 
or blockage and the zero setting of the transducer remained constant the 
dynamometer was started again, brought to the required test speed and a pressure 
scan was made. Following completion of the scan, the dynamometer was stopped 
and the mylar tape was removed from the next station to be tested and the last 
tested station was blocked with tape. The dynamometer was then started again 
and another pressure scan was made. Testing continued in this manner until all 
spanwise stations on each blade had been measured. 

Data Reduction 

Pressure measurements were recorded by the pressure transducer at the 
center of rotation of the hub assembly. Hence a correction must be applied to 
account for centrifugal effects along the brass tubing leading to the pressure 
search holes. The pressure gradient is given by 


SL where V - 


From the characteristic gas equation p 55 pRT 
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Substituting equation (9) into equation (8) and integrating 


locu p - / 4- K. 

J 1 (^£ 7 , 


( 10 ) 


When r=0, p-p^measured valve of the pressure s hence K 


Hence 


F ; P« 


f-li V' 

,vyv r 


log p 
e e m 


(ID 


(u^angular velocity in rads/sec; R=the characteristic gas constant in ft.lb/lb/°F 

2 

T®temperature in°R; g=gravitational constant, 32,2 ft/sec ; r=radius in feet) 

This pressure is perpendicular to the surface of the propellor. Also the 
measured pressure must be corrected by any initial pressure measured by scanivalve 
tubulations located along the axis of spin in the hub. Hence the pressure 


perpendicular to the axis of rotation is 


?= (f~- 


e 


a 


cr£ 


where p - pressure correction and © is the local chord angle * The transducer 

- - i>,- 

calibration was determined to be 0* 001225 volts per square foot of pressure. 

■j. 

All results were measured with the propellor rotating at a tip speed of 500 ft/sec 


RESULTS 

Typical pressure distributions at chord positions given by r/R = 0-3,0.5,0.75 
and 0.95 are shown in figures 6 and 7. Spanwise net load distributions are shown 

- 5, lO^and 15? This figure clearly demonstrates that 
the load distribution is a maximum near the tip radius and increases with the 
angle. Integration of figure 8 would give the total load. 
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Data Acquisition Block Diagram 
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